Marginality of bulk-edge correspondence for single-valley Hamiltonians 
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We study the correspondence between the non-trivial topological properties associated with the 
individual valleys of gapped bilayer graphene (BLG), as a prototypical multi- valley system, and the 
gapless modes at its edges and other interfaces. We find that the exact connection between the 
valley-specific Hall conductivity and the number of gapless edge modes does not hold in general, 
but is dependent on the boundary conditions, even in the absence of intervalley coupling. This non- 
universality is attributed to the absence of a well-defined topological invariant within a given valley 
of BLG; yet, a more general topological invariant may be defined in certain cases, which explains 
the distinction between the BLG-vacuum and BLG-BLG interfaces. 



I. INTRODUCTION 

The analysis of the electronic properties of semicon- 
ductors often requires considering the presence in their 
band structure of multiple degenerate valleys, centered 
around different, symmetry-related positions in recipro- 
cal spaced At low energy, the dynamics of electrons in 
each of these valleys can be modeled in terms of a long- 
wavelength effective Hamiltonian, an approach suitable 
to describe the vast majority of the transport phenom- 
ena that are measured experimentally. One interesting 
question - which is particularly timely in view of the 
intense research effort focusing on so-called topological 
insulators^! - j s whether there exist non-trivial topolog- 
ical properties associated to the individual valleys that 
can be described in the framework of a long-wavelength 
effective Hamiltonian. 

To start addressing this question, here we focus on the 
case of two-dimensional (2D) electronic systems. The 
two best known examples of non-trivial topological insu- 
lators in two dimensions are provided by integer quan- 
tum Hall sys tems^ and by spin-orbit induced topological 
insulators^'. In these systems, the non-trivial topologi- 
cal properties of the bulk band structure result in a quan- 
tized Hall conductivity when the Fermi level is located 
in a bulk energy gap (for spin-orbit induced topological 
insulators, spin-Hall conductivity is quantized only when 
the spin is a good quantum number and one can consider 
the Hall conductivity for each spin state separately). In- 
deed, it is well established that the expression of the bulk 
Hall conductivity an given by the Kubo formula for lin- 
ear response corresponds to the Chern number that char- 
acterizes the topological structure of the mappin g fr om 
the Brillouin zone to the space of the Bloch states-^. If 
we confine ourselves to the simplest case where only two 
bands are relevant, the physics of these systems can be 
described in terms of an effective Hamiltonian of the form 



H = 



k 



#[[-g(fc) • a + e k -fj]^ k 



(1) 



where — [ck^,Cki\ T is the itinerant electron field op- 
erator, <j = (cr x , a y , <t z ) is a vector of the Pauli matrices, 




FIG. 1: Plot of B k = • [9 fex g x d kv g] (integrand of Eq. 
(J3|) as a function of momentum k over the entire Brillouin 
zone, obtained from the tight-binding description of gapped 
BLG. Bk is nonvanishing only close to the K and K' points 
(valleys), where it is well approximated by the expression for 
g(fc) that enters the long-wavelength effective Hamiltonian 
(see Eq. Q). The integral of B h over the entire Brillouin 
zone vanishes, owing to the equal and opposite contributions 
of the two valleys (= ±1, corresponding to the quantized Hall 
conductivity of the individual valleys). 



and g(fc) is a real vector. If the bulk band structure is 
fully gapped, and the chemical potential fj, lies inside the 
gap, the Hall conductivity (in units of e 2 /h neglecting 
spin degeneracy), equal to the Chern number associated 
to the occupied band, is given by^ 



N 



^ J d 2 kg- (<9 fex g x d ky g) 



(2) 



with g = g/|g| (again, for spin-orbit induced topologi- 
cal insulators with conserved spin, these expressions hold 
separately for the two spin directions). The integral has 
to be performed over the entire Brillouin zone. The sys- 
tem is topologically non-trivial if N ^ 0. The physical 
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manifestation of the topological non-triviality, which dis- 
criminates between topological and trivial insulators, is 
the appearance of N gapless states at the system edges, 
which can transport current even when the Fermi level is 
located in the bulk energy gap. This statement, true for 
both integer quantum Hall and spin-orbit induced topo- 
logical insulators with conserved spin, relates the bulk 
electronic structure (which determines the Chern num- 
ber) to the e dge p roperties and is known as bulk-edge 
correspondence^^. 

The integrand in Eq. ^ is often sharply peaked in the 
regions of the momentum space where |g(k)|, and hence 
the gap, becomes small - these are the valleys of the band 
structure. Thanks to the fast convergence, the integral 
m Eq. (§ over the Brillouin zone can be calculated using 
the long wavelength approximation for the Hamiltonian 
valid in the individual valleys. In this case one can as- 
sociate an index N T to each one of the valleys (r labels 
the valley), and the Chern number is then the sum of the 
N T . This sum can vanish even if the individual N T do 
not, in which case the system is -according to the def- 
inition given above- topologically trivial. However, the 
question arises as to whether the non- vanishing of the N T 
associated to the individual valleys has observable conse- 
quences. Indeed, a non-vanishing N T implies that valley 
r gives an non-zero contribution to the Hall conductivity, 
and the appearance of N T gapless edge states associated 
to each individual valley may be expected by a naive 
"extension" of the bulk-edge correspondence. It may be 
argued that due to the overall triviality of the system 
(i.e., the fact that the Chern number defined over the en- 
tire Brillouin zone vanishes) these states localize because 
they are not protected against disorder at the edge, which 
causes inter- valley scattering^. Nevertheless, the funda- 
mental question remains, whether in a system with ideal 
edges (i.e., edges that preserve the valley quantum num- 
ber, well defined at low energies) there exists a bulk-edge 
correspondence for individual valleys. 

Here we consider the case of gapped bilayer graphcnc 
(BLG) as a model system for the case in which two 
symmetry-related valleys are present (the K and K' val- 
leys with non- vanishing iV r = ±1 depending on the val- 
ley, see Fig. [I]), and investigate the low-energy electronic 
states at its edges. Specifically, we consider different 
crystalline edges which do not couple the valleys, and 
we find the corresponding edge states in the different 
cases (we have discussed elsewhere the effect of disor- 
dered edges, which is experimentally more relevant!^). 
If the bulk-edge correspondence could be generically ex- 
tended to single-valley Hamiltonians, one would expect 
that N T — ±1 should imply the presence of exactly one 
gapless mode per valley per spin at an ideal edge of 
gapped BLG, for all edges that do not couple the valleys. 
In contrast to this expectation, we show through an ex- 
plicit analytical solution in complete agreement with full 
tight-binding calculations, that the number of the gapless 
edge modes depends on the boundary conditions, even 
when no mixing between valleys exists at the edge^. In 
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FIG. 2: Illustration of bilayer graphene edges (highlighted by 
thick lines) terminating all in the zigzag direction but with 
different sublattices. The number of subgap edge modes is 1 
for case (a) and (b), 2 for case (c) and for case (d). A unit 
cell of the lattice is also shown in (a) in broken lines. 



other words, the valley- specific bulk-edge correspondence 
is not fulfilled for plain BLG-vacuum interfaces. 

In order to understand this result, we examine the ge- 
ometrical meaning of N T associated to a particular valley 
Hamiltonian, and find that, contrary to the case of the 
Chern numbers in topological insulators (defined over the 
entire Brillouin zone), N T in BLG does not correspond in 
itself to a well-defined topological invariant of a mapping. 
Nevertheless, the non-vanishing of N T still signals that 
the properties of the electron states in individual valleys 
are non-trivial, and we discuss how this non-triviality 
manifests itself at interfaces between different domains 
where the gap of BLG changes sig rP. In this case, the 
number of zero-energy states corresponds to the differ- 
ence of iV T on opposite sides of the BLG-BLG interface 
(iV* — N£, where we denote with N l T ' r the value of N T at 
the left and at the right of the BLG-BLG interface), in 
agreement with expectations based on bulk-edge corre- 
spondence. Indeed, we show that for a BLG-BLG inter- 
face the difference N l T — is a well-defined topological 
invariant of a mapping, even though the individual N l T 
and N£ are not. These results reveal the marginal char- 
acter of individual valleys in gapped BLG: the difference 
of N T across an interface may or may not be topologically 
well defined, depending on the specific kind of interface. 
It is a well defined topological invariant across a domain 
wall where the gap changes sign, but it is not for a BLG- 
vacuum interface. Such a marginal topological character 
of a single valley is likely to be not only characteristic of 
BLG, but is a more general property common to many 
other multi-valley systems. 
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II. ABSENCE OF VALLEY-SPECIFIC 
BULK-EDGE CORRESPONDENCE 

We start by briefly reviewing the known properties of 
the band structure of bilayer graphene, emphasizing the 
specific aspects that will be relevant later. In BLG a gap 
between valence and conductance bands can be opened 
controllably by applying a perpendicular electric field, 
while maintaining the Fermi energy in the middle of the 
g ar JHHl2l (th e s jg n f th e g a p j s determined by the ori- 
entation of the field). The continuum low-energy, long- 
wavelength limit Hamiltonian for each individual valley 
can be obtained from the tight-binding description of 
BLG with Bernal stacking^. With the four inequivalent 
atoms in a unit cell labeled by A(B)1(2) (see Fig. [2ja)), 
the dimcnsionlcss Hamiltonian matrix for the wave func- 
tion * = (xBi,XA2,fB2,<PAi) T readsPl 

Ht = ( k\ Ho)' (3) 
H i = a x + Aa z , H = -Aa z 
K T = K\ = rk x a x + k y a y . 

Here, k = (k x ,k y ) is the (dimensionless) wave vector, 
the valley index t = ±1, and 2|A| (|A| <C 1) defines the 
size of the bulk gap. The excitation spectrum is E^i ~ 

±^A 2 + (k 2 + ( 2 ll ) 2 , witn eh = i an d £; = o. This 

implies a splitting of an order of magnitude 1 between 
the high energy bands (h) and direct gap 2|A| between 
the low energy bands (I). 

The original four-component Hamiltonian (|3| can be 
inconvenient to work with, especially when one needs to 
find solutions with specific boundary conditions. There- 
fore, bearing in mind that the energy range of our in- 
terests would only involve the two low energy bands, we 
reduce (J3| to a two-component Hamiltonian by rewriting 
the Schrodinger equation H T ^ = in the following 
form: 

X = -(E + a x + Aa z )K T <p, (4) 
-\K T a x K T + A(l - k 2 )a z ]ip = E(l + k 2 )ip, (5) 

where <p = ((fiB2,<PAi) T and x = (xbi,XA2) T - As we 
are interested in the subgap edge state solutions with 
\E\ < |A| <C 1, only terms up to linear order in E and 
A are kept. If we further neglect 0(Ak 2 ) and OiEk 2 ) 
terms, we obtain 

/ A (rk x — ik v ) 2 \ „ 

-{(Tk x+i k v r -a v> )*= E ^ ^ 

which is the reduced Hamiltonian, describing the low en- 
ergy bands to the lowest order in A and k. This re- 
duction procedure is essentially identical to that used by 
McCann and Fal'kcP. Here we emphasize two important 
aspects: first, the nontrivial topological properties of the 
original four-component model are fully inherited by the 
reduced two-component model - this will be shown in 



Section |H A| second, the relation between the low-energy 
(ip) and high energy (x) components of the wavefunc- 
tion [Eq. Q] is essential to impose the correct boundary 
conditions at edges - this will be shown in Section |II B| 

A. Quantized Hall conductivity of individual 
valleys 

Starting from the effective Hamiltonians given previ- 
ously and following the logic explained in the introduc- 
tion, we examine the nontrivial topological properties of 
a single valley of gapped BLG by calculating the corre- 
sponding (single-valley) contribution to the Hall conduc- 
tivity <r T H . We consider the case when the Fermi energy 
is lying in the middle of the bulk energy gap, so that we 
are dealing with one completely full and one completely 
empty bulk bands of The valley-specific Hall conduc- 
tivity ajj (in units of e 2 /h) corresponds, through Kubo 
formula, to the quantity N T defined previously. Here 
we present calculations for both the original four-band 
model ([3| and the reduced two-band model ^ , and show 
that these calculations yield consistent results. This con- 
sistency implies that the nontrivial momentum topology 
of the original four-band model is fully inherited by the 
two-band model, which lends itself more easily to a geo- 
metrical interpretation of our results. 

The Hall conductivity (in units of e 2 /h) assoc iated to 
a fully-gapped model is generically given b} * 20 * 21 ! 

N= iS Tr / < P kdk o G K G ~ lGd *u G ~ lGd kiG- 1 , 

(7) 

where the non-singular propagator G(k ,k) is defined as 
G = (iko — H)^ 1 . For the four-component single- valley 
Hamiltonian of BLG p]), we find iV r = rsign(A), which 
is equal and magnitude but opposite in sign (±1) for two 
valleys; therefore, the Chern number defined over the 
entire Brillouin zone vanishes. This result is expected 
for BLG, which is strictly speaking topologically trivial, 
despite the fact that the individual valleys possess non- 
trivial N T . 

For the reduced two-component model given by Eq. 
([Tj]), we can write the Hamiltonian as H T 2x2 ^ = — g r (fc) • 
<r, where g r (fe) = {k 2 — k 2 , 2rk x k y , A). In this case, 
using the simplified form of the above formula for two- 
band models^, given by ([2|, we find N T = rsign(A), 
in agreement with the calculation based on the full four- 
band model. This agreement indicates that the nontrivial 
topological properties of the individual valleys are due 
solely to the low energy bands. 

It is worth noting that the integral (divided by 2) in 
Eq. (§ can be also identified as the Berry phase acquired 
by an electron adiabatically transferred along the path 
enclosing the area of the integraP^. In this sense, the 
quantized Hall conductivity is equivalent to the quantized 
Berry phase (in units of 27r) associated with the single- 
valley HamiltoniarP^. 
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FIG. 3: Comparison of the dispersion relations for the subgap 
edge modes obtained by solving the continuum model and 
exact diagonalization of the tight-binding model, with A = 
0.15. The upper panel case corresponds to the edge shown 
in Fig. [5] (a), and contains only one branch of gapless edge 
modes; the lower panel case corresponds to Fig. [2] (c), and 
contains two branches of gapless edge modes. 



B. Boundary conditions and zero-energy edge 
states 



We now proceed to calculate the exact solution of the 
edge states in gapped BLG, for different edge structures 
that do not mix the valleys. We will only consider one 
of the two valleys, with r = +1, and the results for the 
other valley can be inferred through a symmetry trans- 
formation similar to time reversal. 

The generic constraint for boundary conditions in the 
present model is imposed by vanishing probability cur- 
rent across the boundary. This constraint allows for a 
variety of boundary conditions which are physically valid 
(see Appendix [V]) . In order to base our investigation on 
concrete physical examples, we illustrate the boundary 
conditions that will be discussed in the following in the 
tight-binding picture of BLG. Fig. [2] shows the edge 
structure for four cases, where the BLG lattice termi- 
nates at different sublattices, either A or B, in each of 
the two layers. Indeed, each of these combinations corre- 
sponds to a distinct boundary condition imposed in the 
continuum model, similar in spirit to what has been dis- 
cussed by Brey and Fertig^lfor single layer graphene. A 



significant difference is that in BLG the boundary condi- 
tions in general involve both low (ip) and high (x) energy 
components of the wavefunction. Still, thanks to Eq. Q, 
the latter can be written in terms of the low energy com- 
ponents and their derivatives, which is why a description 
in terms of the low energy components only is possible. 
For the structure in Fig. [2ja) where the BLG lattice 
terminates at the edge with sublattices Al and A2, this 
corresponds to imposing tpB2(y = 0) = XBi(y = 0) = 
(the BLG is infinitely long in the x direction and semi- 
infinite in the +y direction, i.e. the edge is located at 
y = 0). The general solution of §5§ for \E\ < |A| is 



with 



K 2 



E 



(A+E)- 



1 

-(A- 



E)(kl-n 2 B ) 



ik x w—K s y 



(8) 



0*+K S ) 2 



k\ - (A 2 + E 2 ) ± i VA 2 - E 2 , »(«±) > 0, (9) 



and c s (s — ±) are coefficients that need to be deter- 
mined. Using Q, the boundary conditions lead to two 
coupled equations for c s with the secular equation for the 
existence of nontrivial subgap solutions 



{k x + K + )(k x + «_) = {A + Ef 



(10) 



This equation gives the dispersion relation E(k x ) for one 
sub-gap edge mode. When \E\, |A| <C k 2 (as is true for 
the E — mode) Eq. |To| ) can be simplified (since n± ~ 
\k x \±i(VA 2 ~E 2 )/2\k x \) to: 



E 
A 



1 + m ' 



k x < 0. 



(ii) 



In Fig. [3](a) the analytical solution is compared with 
the full tight-binding solution, from which we see that 
the continuum model represents an excellent approxima- 
tion. We conclude that for the edge just considered, one 
zero-mode exists (per valley and spin^PU, i.e. the result 
expected if the bulk-edge correspondence for single valley 
holds. 

We now consider the edge shown in Fig. [2jc). The 
corresponding boundary conditions are <PB2\JJ = 0) = 
<PAi{y — 0) = 0, which lead to the new dispersion relation 
for the subgap edge modes, given by 



*2k x {k x 



-)(** 



A + E 
A-E' 



(12) 



This equation has two solutions related by transforma- 
tion k x — > — k x and E — > —E, which comes from the 
symmetry of the wave equation (|5| under {k x , E, tp} — > 
{~k x ,—E,ia y ip}. It is this symmetry, which is also pre- 
served by the boundary conditions in the present case 
(whereas it is broken by the boundary conditions con- 
sidered in the previous example), that guarantees the 
existence of two subgap edge modes. The dispersion 
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relation given by (12 1 can be further simplified when 
|£|,|A| < hi, to 



A=~ sgn(fca;) TT^!' 



(13) 



which explicitly shows a pair of gapless edge modes prop- 
agating in the same direction, plotted in Fig. [3jb) 
(again, in excellent agreement with the tight-binding so- 
lution). Without presenting the details of the solutions 
for the other cases, we only state that the edge shown 
in Fig. §» (<p A1 (y = 0) = X A2(y = 0) = 0) leads 
to a single branch of subgap edge modes, similar to the 
first example above, and the one shown in Fig. [2jci) 
{XBi(y = 0) = XA2(y = 0) = 0) yields no subgap modes. 
We therefore have to conclude that the number of zero- 
modes at the BLG- vacuum interface depends on the spe- 
cific edge considered, i.e., the bulk-edge correspondence 
relating the Hall conductivity and gapless edge modes 
does not hold in general for individual valleys^. 



III. DISCUSSION 

Finding that bulk-edge correspondence does not hold 
for an individual valley, even though the contribution 
that each valley gives to the Hall conductivity is inte- 
ger (in units of e 2 /h) may seem surprising. Specifically, 
one may wonder why the argument based on Laughlin's 
gedanken experiment!^ - which is normally invoked to 
justify the existence of edge states in integer quantum 
Hall systems - does not apply to individual valleys in 
gapped BLG. The answer to this question has to do with 
the fact that Laughlin's gedanken experiment considers 
the adiabatic flow of electric charge, which is a conserved 
quantity (in the context of Laughlin's argument conser- 
vation of charge is so obvious that this assumption is not 
normally emphasized). Contrary to the charge, the valley 
quantum number is not in general conserved. In partic- 
ular, it is not conserved even in the presence of ideal 
edges if one considers the process - the adiabatic inser- 
tion of flux, which is the basis of Laughlin's gedanken 
experiment. As we outline below, the non-conservation 
of the valley quantum number allows for the possibility 
of a quantized valley Hall conductivity in the absence of 
zero-energy edge modes. 

Upon varying the magnetic flux (j> in the usual cylin- 
drical geometry involved in Laughlin's argument, the al- 
lowed k x values of momentum along the edge flow ac- 
cording to k x — > k x + 4>/L, where L is the circumference 
of the cylinder. Upon insertion of flux </> = 2ir, the set of 
allowed k x values returns to the original one. From the 
generic form of the electron wave function Q it follows 
that the wave-function weight is redistributed in the di- 
rection perpendicular to the edge in the process of flux 
insertion. The redistribution occurs in the opposite di- 
rections for the states near two valleys, with states in one 
valley flowing towards the edge, and states in the other 
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FIG. 4: Illustration of the mapping g(fc) for gapped BLG. 
Panel (a) shows g(fc) for a single valley as a vector field on the 
R 2 plane of k. Panel (b) shows the solid angle Q that g covers 
when k runs over the whole plane, the double arrow signifying 
the double covering of the (upper) hemisphere. Panel (c) 
shows how two marginal topological mappings can be "glued" 
to form a well-defined topological mapping, which is the case 
for a BLG-BLG domain wall with opposite mass signs on the 
two sides. 



valley flowing away from the edge. This is the micro- 
scopic mechanism for the transverse valley current flow 
implied by the non-zero value of the valley Hall conduc- 
tivity, a v H = N T — N- T . The fact that the valley current 
flows perpendicular to the edge, according to continuity 
equation would imply accumulation of the valley density 
at the edge. This apparent valley charge accumulation 
can be accommodated in two ways. The first is by creat- 
ing valley charge imbalance at the edge, in direct analogy 
to quantum Hall and quantum spin Hall effect with con- 
served spin. Clearly, for that to occur, the gapless edge 
modes have to exist. Alternatively, the valley current 
influx at the edge can be compensated by the adiabatic 
transfer of the electrons between the valleys induced by 
the flux insertion, which converts k x — ¥ k x + 2n/L for 
states in the entire ID Brillouin zone of the cylinder. 
This is possible due to the fact that the valleys are in- 
deed connected by the bands that extend deep below the 
chemical potential. As a consequence, the spectral flow 
cannot be fully described in terms of the long- wavelength 
effective Hamiltonian in the neighborhood of the K and 
K' points, but it requires considering the Hamiltonian of 
the system over the entire Brillouin zone. This is how the 
non-conservation of the valley quantum number can ac- 
count for the seemingly paradoxical result that non-zero 
a v H can exist even when gapless edge modes are com- 
pletely absenlpSl. 

Next, let us give a geometrical reason for the the ab- 
sence of a valley-specific bulk-edge correspondence at 
gapped BLG edges (BLG-vacuum interfaces). It origi- 
nates from the fact that the quantity N T associated with 
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a given valley is not a well-defined topological invariant. 
This statement may appear to be in conflict with the 
case of domain wall interfaces, across which the sign of 
the mass (A) in BLG changes sigrP^. In this latter case, 
bulk-edge correspondence "predicts" that gapless states 
should be present, whose number corresponds to iVj— iV£. 
When the gap A changes sign across the domain wall, 
\N~. — N£\ — 2 and indeed two gapless states are found 
at the interface. The conflict is only apparent, because 
N l T — does correspond to a well-defined topological 
invariant, even though N T does not. 

To understand the difference between N T and N l T — 
we look at the integrand function that is used to calcu- 
late these quantities, which, from Eq. ([2]), is given by 
®k = j^g' [<9fc x g x dk g]- This expression corresponds to 
the Jacobian of the transformation from the region R 2 
of the momentum space in the vicinity of a given val- 
ley point, to the unit sphere S 2 on which g resides. For 
BLG, at large k (k 2 |A|), the vector g lies on the 
equator, and therefore the integral of Bk does not repre- 
sent a topological invariant since the mapping between 
non-compact R 2 and compact S 2 is topologically triv- 
ial. Indeed, as shown in Fig. [4] (a) and (b), the vector 
g wraps around half of the unit sphere twice, but does 
not cover the entire sphere. That is despite the fact that 
the integral rapidly converges for k 2 |A| and is equal 
to 1 (for a given valley and sign of the gap A; it is —1 
for the opposite valley or for opposite sign of A). For 
massive Dirac fermions such non-compact behavior has 
been discussed by Volovik, and is known as marginal. 
In the same sense, the massive chiral fermions in BLG 
are marginal as well. The marginality implies that a 
small perturbation in the Hamiltonian (e.g. momentum- 
dependent mass term A(k 2 )) or boundary conditions can 
have a big effect on presence or absence of the gapless 
modes (e.gP^). 

Now, for a domain wall such that A(x < 0) < and 
A (a; > 0) > 0, despite the marginal character of N T on 
the two sides of the domain wall, their difference is a well- 
defined topological invariant. That is because at k — > oo, 
the A in the Hamiltonians becomes irrelevant and the 
textures of the g vectors for the two insulators seamlessly 
connect on the equator of S 2 , thereby compactifying the 
momentum space within the same valley. The connection 
of the textures is shown schematically in Fig. [4] (c) . Sim- 
ilar considerations are not only valid for BLG, but can 
also be extended to individual valleys in (gapped) single 
layer graphene^, as well as to interfaces between gapped 
single layer graphene and Kanc-Mclc topological insula- 
tors with conserved spirP. In these cases, the differences 
in the single-valley N T across the interface is a well de- 
fined topological invariant, corresponding to number of 
gapless states present (or more precisely, the difference 
between the numbers of left and right moving gapless 
modes in a given valley). Mathematically, this is a con- 
sequence of the index theorem discussed in this type of 
contexts by Volovik (see Section 22.1.4 irP). 



IV. CONCLUSION 

Our work shows that the long wavelength Hamiltoni- 
ans that describe the low-energy electronic states in indi- 
vidual valleys of gapped BLG possess non-trivial topolog- 
ical features. In contrast to integer quantum Hall systems 
and quantum spin Hall systems (with conserved spin), 
these non-trivial valley-specific properties cannot be de- 
scribed by a topological invariant. Rather, they require 
the analysis of the specific interfaces (in the present case 
the BLG- vacuum interface and the domain wall) in order 
to establish whether a topological invariant that deter- 
mines the number of gapless edge modes can be defined. 
This is the characteristic signature of marginal topologi- 
cal insulators. 

It is certainly the case that the low-energy edge states 
that can be predicted through these topological consid- 
erations are not robust against short-range disorder that 
couples the valleys (even though they are stable against 
long range disorder that scatters electrons within the 
same valley). In the presence of inter- valley scattering, 
the gapless modes associated to individual valleys will 
couple, leading to the opening of a transport gap. Nev- 
ertheless, this does not mean that the low-energy states 
at interfaces and domain walls are experimentally irrel- 
evant. In fact, for sufficiently pure materials when the 
Fermi energy is in the gap, the states originating from 
the localization of the gapless modes can provide a dom- 
inant path for transport in the insulating state, because 
no other states are available deep in the bulk gap of the 
material. 

Finally, even though all calculations presented here 
have been performed specifically for gapped BLG, our 
arguments do not rely on the particular form of the val- 
ley Hamiltonian, or the number of valleys. Consequently, 
similar treatment based on other effective valley Hamil- 
tonians should be applicable to many other materials and 
interfaces, and may provide a convenient tool for deter- 
mining the interfacial electronic properties. 
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V. APPENDIX: BOUNDARY CONDITIONS 

In this appendix we derive the most general boundary 
conditions from the continuum model given by pi or 
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([6]), by imposing current conservation through a hard- 
wall boundary. We show that the boundary conditions 
associated to the edge structures considered in the main 
text are specific realization of the general case. 

Let us start from the original four-component model 
([3]), for which the current operator is given by 



3 = 



a 
a 



(14) 



where er = (o Xl o y ) and we have let h = 1. Without 
losing generality, we consider a boundary at y = with 
unit normal ris = — y; then, the current conservation 
implies 



i(Xm<PAi 



X*A2 ( PB2)\ y=0 + C.C. 



= 0. 

(15) 



We immediately see that the above condition can be nat- 
urally satisfied by letting in each layer the wave ampli- 
tude on either A or B component vanish at y = 0. This 
leads to four choices of combinations each finding an ex- 
act correspondence to the cases listed in Fig. [2] Indeed, 
the same amount of freedom in "choosing" boundary con- 
ditions is also inherited by the reduced model ([6|, for 



which, by the same token, the vanishing current across 
the boundary requires 

{ip* A1 d+ip B 2 - <p B2 d-(p A i)\ y=0 + c.c. = 0, (16) 

where d± = d x ± id y . Noticing that the leading order 
approximation of Eq. Q yields 



XBl = id + (fB2, XA2 = id-ipAl, 



(17) 



Eq. (16) is nothing but Eq. (15). In other words, the 



reduced model shares the same choices of boundary con- 
ditions as the original one, and the edge modes it allows 
for may vary accordingly. 

These derivations indicate that for a generic multi- 
component Hamiltonian, considerable freedom exist in 
selecting boundary conditions associated to the wave- 
equation, reflecting different possible structures of the 
edges. To demonstrate the dependence of the gapless 
edge modes on boundary conditions, which serves the 
purpose of the current work, however, the four cases pre- 
sented in the main text are already sufficient. 
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